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Abstract: In two recent works O [T], it has been shown that the counting generating 
functions (CGF) for the 23 walks with small steps confined in a quarter-plane and 
associated with a finite group of birational transformations are holonomic, and even 
algebraic in 4 cases - in particular for the so-called Gessel's walk. It turns out that 
the type of functional equations satisfied by these CGF appeared in a probabilistic context 
almost 40 years ago. Then a method of resolution was proposed in j^, involving at once 
algebraic tools and a reduction to boundary value problems. Recently this method has 
been developed in a combinatorics framework in [TT] , where a thorough study of the explicit 
expressions for the CGF is proposed. The aim of this paper is to derive the nature of the 
bivariate CGF by a direct use of some general theorems given in 
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Sur I'holonomie ou I'algebricite de fonctions generatrices 
comptant des marches aleatoires dans le quart de plan 



Resume : Dans deux articles recents [21 [T] , il a ete montre que les fonctions generatrices 
de comptage (CGF) pour 23 marches avec pas unite, evoluant dans le quart de plan et 
associees a un groupe fini de transformations bi-rationnelles, sont holonomes et meme 
algebriques dans 4 cas - notamment celui de la marche dite de Gessel. Or il s'avere que le 
type d'equations fonctionnelles verifiees par ces CGF est apparu il y a environ 40 ans, dans 
le contexte probabiliste des marches aleatoires. Une methode de resolution fut proposee 
dans ID, basee a la fois sur des techniques algebriques et sur la reduction a des problemes 
aux limites. Recemment cette methode a ete developpee dans un cadre combinatoire 
dans ou une etude fouillee des formes explicites des CGF a ete menee. L'objet du 
present article est de retrouver la nature des CGF bivariees par une application directe de 
theoremes donnes dans j^. 

Mots-cles : Algebrique, holonome, fonction generatrice, marche discrete homogene par 
morceaux, quart de plan, couverture universelle, fonction elliptique de Weierstrass. 
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1 Introduction 



The enumeration of planar lattice walks is a classical topic in combinatorics. For a given 
set S of allowed jumps (or steps), it is a matter of counting the number of paths starting 
from some point and ending at some arbitrary point in a given time, and possibly restricted 
to some regions of the plane. A first basic and natural question arises: how many such 
paths exist? A second question concerns the nature of the associated counting generating 
functions (CGF): are they rational, algebraic, holonomic (or D- finite, i.e. solution of a 
linear differential equation with polynomial coefficients)? 

For instance, if no restriction is made on the paths, it is well-known that the CGF are 
rational and easy to obtain explicitly. As an other example, if the walks are supposed to 
remain in a half-plane, then the CGF can also be computed and turn out to be algebraic 
in this case, see e.g. [3]. 

Next it is quite natural to consider walks evolving in a domain formed by the intersection 
of two half-planes, for instance the positive quarter-plane Z^. In this case, the problems 
become more intricate and multifarious results appeared: indeed, some walks admit 
of an algebraic generating function, see e.g. [7J and [8j for the walk with step set 
S = {(— 1, 0), (1, 1), (0, — 1)} and starting from (0,0), whereas some others admit a CGF 
which is even not holonomic, see e.g. [3] for the walk with S = {( — 1, 2), (2, — 1)} and 
starting from the point (1,1). 

In this framework, M. Bousquet-Melou and M. Mishna have initiated in [2j a systematic 
study of the nature of the walks confined to Z^, starting from the origin and having steps 
of size 1, which means that S is included in the set {{i,j) : \j\ < 1} \ {(0, 0)}. Examples 
of such walks are shown in figures [TTT] and [T^ below. 

A priori, there are 2*^ such models. In fact, after eliminating trivial cases and models 
equivalent to walks confined to a half-plane, and noting also that some models are obtained 
from others by symmetry, it is shown in [2] that one is left with 79 inherently different 
problems to analyze. 

A common starting point to study these 79 walks relies on the following analytic approach. 
Let q{i,j,k) denote the number of paths in starting from (0,0) and ending at at 
time k. Then the corresponding CGF 



i,j,k>0 



satisfies the functional equation obtained in 



K{x, y)F{x, y, z) = c{x)F{x, 0, z) + c(y)F(0, y, z) + co(x, y), 
K{x,y) = xy 



XI) 



:i.2) 



where 

r 

L {i,i)G5 

c{x)= 

(i ,-l)GS 

c{y)= Yl y'^'' 
(~ij)e5 

cq{x, y) = -(5F(0, 0, z) - xy/z, 

with (5 = 1 if (—1, —1) G 5 [i.e. a south-west jump exists], (5 = otherwise. 

For z = 1/|5|, (|1.2p plainly belongs to the generic class of functional equations (arising in 

the probabilistic context of random walks) studied and solved in the book H], see section 1X1 
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of our appendix. For general values of z, the analysis of (|1.2p for the 79 above-mentioned 
walks has been carried out in where the integrand of the integral representations is 
studied in detail, via a complete characterization of ad hoc conformal gluing functions. 
It turns out that one of the basic tools to solve equation (|1.2p is to consider the group 
exhaustively studied in and originally proposed in [lOj where it was called the group 
of the random walk. Starting with this approach, the authors of |2] consider the group of 
birational transformations leaving invariant the generating function j)eS •'^^y'' 
precisely the group W = {(,,1]) generated by 



X, - 



(« -1)65 



X' 



(j,+i)e5 ■ 



ri{x,y) 



xj2 



(+ij)e5 



• y 



Clearly ^0^ = 7/07^ = id, and is a dihedral group of even order larger than 4. In |^ this 
order is calculated for each of the above-mentioned 79 cases: 56 walks admit an infinite 
group, while the group of the 23 remaining ones is finite. 

It is also proved in p] that among these 23 walks, W has order 4 for 16 walks (the ones 
with a step set having a vertical symmetry), W has order 6 for 5 walks (the 2 at the left 
in figure 11.11 and the 3 ones at the left in figure II. 2p , and W has order 8 for the 2 walks on 
the right in figures 11.11 and 11.21 Moreover, for these 23 walks, the answers to both main 
questions (explicit expression and nature of the CGF (jl.ip ) are known. In particular, the 
following results exist. 

Theorem 1.1 ([2j). For the 16 walks with a group of order 4 md for the 3 walks in 
figure \1.1\ the formal trivariate series is holonomic non-algebraic. For the 3 walks 

on the left in figure the trivariate series U.l\) is algebraic. 

The proof of this theorem relies on skillful algebraic manipulations together with the 
calculation of adequate orbit and half- orbit sums. 

Theorem 1.2 (tlj)- For the so-called Gessel's walk on the right in figure the formal 
trivariate series lil.l]) is algebraic. 

The proof given in [T] involves a powerful computer algebra system (Magma), allowing to 
carry out dense calculations. 

The main goal of our paper will be to present another proof of theorems 11.11 and 11.21 for 
the bivariate generating function (x, y) 1-^ F{x, y, z), by means of the general and powerful 
approach proposed in the book |^. In particular, we are going to show how chapter 4 of 
this book (which deals with walks having a finite group) yields rather directly the nature 
of the bivariate generating functions (jl.ip for the 23 walks associated with a finite group. 
The rest of the paper has a simple organization. Indeed, section [2] recalls the results of H], 
used in Section |3l where we find the nature of the counting generating functions coming 
in (O). 
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Fig. 1.1: On the left, 2 walks with a group of order 6. On the right, 1 walk with a group 
of order 8 
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Fig. 1.2: On the left, 3 walks with a group of order 6. On the right, 1 walk with a group 
of order 8 



2 Some results of jjl] with miscellaneous extensions 

In this section we use the notation of appendix [Aj where is presented the generic functional 
equation (jA.ip studied in fl]. It clearly contains (jl.ip as a particular case, since the 
coefficients q and go occuring in (lA.ip may depend on the two variables x, y. Moreover, 
although this be not crucial, the unknown function ■n{x) [resp. tt (y)] is not necessarily equal 
to 7r(x,0) [resp. 7r(0, y)]. 

Chapter 5 of [4j explains how to solve (lA.ip by reduction to boundary value problems of 
Riemann-Hilbert-Carleman type, and explicit expressions for the unknown functions are 
obtained. 

Independently, when the group of the walk is finite^ chapter 4 of [4j provides another 
approach allowing to characterize completely the solutions of (jA.ip . and also to give 
necessary and sufficient conditions for these solutions to he rational or algebraic. We 
shall proceed along these lines in the sequel. 



2.1 The group and the genus 

Let C (x, y) be the field of rational functions in (x, y) over C. Since Q is assumed to be 
irreducible in the general case, the quotient field C(x,y) denoted by CQ{x,y) is also a 
field. 

Definition 2.1. The group of the random walk is the Galois group T-L = {(,,1]) of 
automorphisms of CQ{x,y) generated by ^ and r] given by 

i(x,y)=\x,-^ • , r][x,y)=\-^ . 

Let 

6 = rj^. 

Then H has a normal cyclic subgroup T-Lq = {6^,i € Z}, which is finite or infinite, and 
Ti/T-Lo is a group of order 2. When the group Ti is finite of order 2n, we have 5" = id. We 
shall write /q, = a{f) for any automorphism a & Ti and any function / € CQ{x,y). 
In [2], the group W is introduced only in terms of two birational transformations. The 
difference between the two approaches is not only of a formal character. In fact W is 
defined in all of C^, whereas in |2] would act only on an algebraic curve of the type 

l{x,y)eC'^:R{x,y) = xy ^ r^jx*?/^ - 1/z =oL 

An immediate question arises: do the notions of order and finiteness coincide in the 
respective approaches? In general, the answer is clearly no. Indeed (adapting [3], lemma 
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4.1.1) for the group to be of order 4 it is necessary and sufficient to have 



0. 



and this condition depends on while W is independent of z. On the other hand, if W is 
finite so is "H, and conversely if % is infinite so is W . In addition, with an obvious notation, 

Order{n) < Order(W). (2.1) 

In the sequel, we shall encounter groups of order not larger than 8. 

Besides, in [3], the algebraic curve is associated with a Riemann surface which is of genus 
g = (the sphere) or g' = 1 (the torus). Accordingly, the universal covering of this surface 
is 

• the Riemann sphere if g = 0; 

• the finite complex plane C^^; if g = 1. 

All 23 random walks we consider in this paper correspond indeed to g = 1, and an efficient 
uniformization is briefly described in appendix [B] Moreover, they have a finite group, say 
of order 2n (this will be made precise in section |3]) . 

2.2 Algebraicity and holonomy 

For any h G CQ{x,y), let the norm N{h) be defined as 

n-l 



N{h)'^l[hs.. (2.2) 
Written on Q{x,y) = 0, equation (jA.ip yields the system 

(2.3) 




where _ _ 

f='^ ^ ^ _ (*))»? 

It turns out that the analysis of (lA.ip given in [3] deeply differs, according N{f) = 1 or 
A^(/) 7^ 1. In section 131 as for the above-mentioned walks, it will be shown that we are in 
the case N{f) = 1 and we need the following important theorem stated in a terse form. 

Theorem 2.2. If N{f) = 1, then the general solution of the fundamental system (E? 
has the form 



w 

■K = Wl+W2^ , 

r 



where 

- the function wi and r are rational; 

- the function W2 is given by 



n-l 
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- on the universal covering Coj, we have 

where (^(w; wi, Wa) =^(^1^3(0;) stands for the classical Weierstrass C,-function (see 19]) with 
quasi periods uji , uj^ ; 

- the function w is algebraic and satisfies the automorphy conditions w = = ws- ■ 

The quantities ui,U2,uj3 have expUcit forms (see [4J, lemmas 3.3.2 and 3.3.3) and the 
automorphism 6 writes simply as the translation 

6{uo) = Lo + uj^, Vo; G Ctj. 

In addition, when the group is finite of order 2n, there exists an integer k relatively prime 
with n satisfying the relation 

nujs = kuj2, 

which in turn allows to prove (see [JJ, lemma 4.3.5) that ^(oj) defined in (|2.5p is not 
algebraic in x{u}). We are now in a position to derive the following corollary. 

Corollary 2.3. Assume the group is finite of order 2n and N{f) = 1. Then the solution 
■n{x) of Hi A . 1\) is holonomic. Moreover it is algebraic if and only if on Q(x,y) = 0, 

n—l I 

E^^ = 0- (2.6) 

fc=o ni=i fs^ 

In fact, the proof is an immediate consequence of the next lemma. 
Lemma 2.4. In the uniformization given in appendix [Bl let us write 

P(^) - 9{x{yj)), 

where g is a fractional linear transform. Then the function W2{p^^{g{x))) defined by i2.4\ ) 
is holonomic in x, where p'^^ denotes the elliptic integral inverse function of p. 

Proof. Since W2 is the product of $ by an algebraic function, it suffices to prove that 
<I>(p~^(g)) is holonomic. It is known that the class of holonomic functions is closed under 
indefinite integration (see e.g. [6j), thus it is enough to prove that [$(p~^(9))]' is holonomic. 
In fact, we are going to show that [$(p~^(g'))]' is algebraic - we recall that any algebraic 
function is holonomic. Indeed, use on the one hand that the derivative of — Ci,3 is Pi, 3, the 
Weierstrass elliptic function with periods a;i,a;3, and on the other hand that 

{p-^\{u) = \l\^u^ - g2U - g^\^l\ 

g2 and g^, being the so-called invariants of p - for these two properties see e.g. [9]. In this 
way we get 

-2nT[^p-\g))]' = — ^- ^[^ipi^3(p-i(<^) -^2/2) +2Ci,3(^i/2)]. 

The function g being rational, g' /[Ag'^ — g2g — g^Y^"^ is algebraic and in order to conclude 
it is enough to prove that pi,3(p~^(5') — (^2/2) is algebraic. Since noj^ = kuJ2, it is shown 
in [U, lemma 4.3.3 that both p and pi.3 are rational functions of the Weierstrass elliptic 
function with periods uji,kuj2- Then it is immediate that pi.3 is also algebraic function of 
p. Moreover, by the well-known addition theorem for the Weiertrass elliptic functions (see 
e.g. O), p is a rational function of p( • -|- W2/2), so that pi^3 is also an algebraic function 
of p( • -1-^2/2). In particular pi,3(p~^(5') —i02/2) is algebraic and lemma fl^ is proved. ■ 
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Of course similar results can be written for n. In particular, it is easy to check that if 
N{f) = 1 then also N{f) = 1, where N{f) = 1 is defined in terms of d = ^77. Thus we are 
in a position to state the following general result. 

Theorem 2.5. Assume the group is finite of order 2n and N{f) = 1. Then the bivariate 
series 7r(x, y) solution of iA.l\) is holonomic. Furthermore 7r(x, y) is algebraic if and only 
if i2. 6\) and \2.0\) hold on Q{x, y) = 0, l\2.6\) denoting the condition for 7r{y) to be algebraic, 
obtained by symmetry from [27l 



3 Nature of the counting generating functions 

We return now to equation (|1.2p and make use of the machinery of section |2] with 

f='A v. = ^-Ml. (3.1) 

Beforehand it is important to note that the variable z plays here the role of a parameter. 
Moreover, the CGF (jl.ip is well defined and analytic in a region containing the domain 

{(x,y,z) gC^ : |x| < l,\y\ < l,\z\ < 1/\S\}. 

The analysis in jl] aimed at finding the stationay measure of the random walk, in which 
case z = 1. But it is possible to check mutatis mutandis that most of the topological 
results remains valid at least for ^(z) > [see for instance [5j or pT] for time dependent 
problems] . 

Lemma 3.1. For any finite n and f given by \3. we have N{f) = 1. 

Proof. Since c [resp. c] is a function solely of x [resp. y], we have c = and c = c^. Hence 

^ _ c _ eg _ 

Cri Cs-q C 

which yields immediately N{f) = 1 from the definition (|2.2p . ■ 

In order to obtain theorems ll.ll and ll.2l for the bivariate function {x,y) 1— t- F{x,y,z), it is 
thus enough, by theorem 12.51 and lemma 13. 1[ to prove the following final proposition. 

Proposition 3.2. For the 4 walks in figure [LE. h2. 6]) and h2. 6)) hold in €? . As for the 19 

other walks (3 in figure [7771 and the 16 ones with a vertically symmetric step set S), \2. b]) 
or i2.6\} do not hold on {{x,y) £ : K{x,y) = j)e5^*y"' ~ 1/^] =0}? '^''^V ^■ 

Proof. The proof proceeds in three stages. 



(!) The walks in figure{T7^ 

Let us check (|2.6p for the 4 walks in figure 11.21 We begin with the popular Gessel's 
walk, i.e. the rightest one in figure 11.21 Here c = 1 in (|1.2p . so that by (|3.ip / = 1 and 
ip = [{xy)rf — xy]/z. Moreover, the order of the group W was shown in [2] to be equal to 8. 
According to (|2.ip . we know that OrderiTi) < 8. On the other hand, by lemmas 4.1.1 and 
4.1.2 of 15, a simple algebra (details are omitted) can show that the order of Ti cannot be 
equal to 4 neither to 6, which entails Order{T-L) = 8, hence n = 4. 
Then 

"■"^ lh 1 1 

=^'^5" = -Yl [(^y\sf^ ~ (^y)s''] =--Yl sign(/9)(xy)p. (3.2) 

A;=0 lli=l •'<5» k=0 k=0 pen 
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The last quantity in (|3.2p can be understood as ( — 1/z) times the orbit sum of xy through 
the group 7i, and it was remarked in [2j that this orbit sum equals zero for Gessel's walk, 
which is equivalent to (|2.6p in C^. 

For the second walk in figure [L2l Order{W) = Order[l-i) = 6 (n = 3) and still c = 1, so 
that the above argument applies. 

For the walk on the left in figure II. 2| we have 

i = {x,l/{xy)), r] = {l/{xy),y), 

hence ^rj = {l/{xy),x). Moreover c = x, thus by using (|3.ip we get / = x'^y and 
i/j = [y — {xy)'^]/z. Next one easily computes fs = l/(xy^), f^'i = y/x, ifjs = [x — l/y'^]/z 
and ■0^2 = \l/{xy) — l/x^j/z. Finally an immediate calculation shows that (|2.6p holds in 

As for the last walk in figure 11.21 we could check n = 3 and again (|2.6p holds in . 

Concerning equation (|2.6p . only Gessel's walk needs some special care, since the 3 others 
give rise to symmetric conditions. In fact, from a direct calculation along the same lines 
as above, it is not difficult to see that (|2.6p holds in C^. 



(ii) The 16 walks with a vertical symmetry. 















A 






















N0 



Fig. 3.3: Three examples of walks having a group of order 4, with symmetrical West and 
East jumps 



We will show the equality 



X 



n-1 



{x - Xn){y - y^), 



(3.3) 



which clearly implies that (|2.6p will not be satisfied on |(x, y) € : K{x, y) = 0|, for any 
z. 

For these 16 walks we have Order(W) = 4 = Order{T-L), that is n = 2. Moreover, the 
vertical symmetry of the jump set S has two consequences : ffist, necessarily ry = (1/x, y); 
second, the coefficient c in (|1.2p must be equal to x, 1 + or 1 + j; + x^. In addition, for 
any of these three possible values of c, we have / = c/c,, = x^. Hence finally 



n-1 



ZCn 



zc{x'^)s 



jjjk ^^j^!h = (^^)»7 - (a:y)g - {xy)s ^ y{xr^ - x) ^ y^{x - Xr^) 

Factorizing by x^/c and using again c/cq = x^, we obtain (|3.3p . 



(iii) The 3 walks in fiaure [7771 

We are going to show that (|2.6p does not hold on K{x,y) = 0, for any z. 
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By similar calculations as in (|3.3p , we obtain that for the 2 first walks in figure 11.11 (which 
are such that n = 3), 

_ {x - y^){l - xy){y - x^) 



fc=0 



nil fs^ ^y't 



with t = y for the walk at the left and t = x + y for the walk in the middle, so that, clearly, 
does not hold on G : K{x,y) = 0|, for any z. 



At last, for the walk at the right in figure [T7T| we have n = 4, and we get 

ipSk (y-l)(x2-l)(x2-y)(x2-y2) 



E 



k=0 

which is obviously not identically to zero in the set G : K{x,y) = O}, for any z. 



A The basic functional equations (see [4], chapters 2 and 5) 

In a probabilistic framework, one considers a piecewise homogeneous random walk with 
sample paths in Z^. In the interior of Z^, the size of the jumps is 1. On the other hand, 
no assumption is made about the boundedness of the upward jumps on the axes, neither at 
(0, 0). In addition, the downward jumps on the x [resp. y] axis are bounded by L [resp. M], 
where L and M are arbitrary finite integers. Then the invariant measure {TTij,i,j > 0} 
does satisfy the fundamental functional equation 



Q{x, y)n{x, y) = q{x, y)iT{x) + g(x, y)-K{y) + ^0(2;, y), 



(A.l) 



with 



'^{x,y) = ^ T^ijx' ^y^ \ 



TT 



(x) = ^TTiox* \ 7r(y) = ^vToj-y^ S 



Qix,y) = xy 
q{x,y) = X 

q{x, y) = y 



1 - ^ Pijx'^y^ 
«>-i,i>o 



'>o,i>-i 
So{x,y) = TTQo[PoQ{x,y) - 1]. 



In equation (lA.ip . the functions 7r(x, y), 7r(x), 7r(y) are unknown and sought to be analytic 
in the region G C'^ : |x| < 1, \y\ < 1} and continuous on their respective boundaries, 

while the q, q, qQ (up to a constant) are given probability generating functions supposed to 
have suitable analytic continuations (as a rule, they are polynomials when the jumps are 
bounded). The function Q{x,y) is often referred to as the kernel of (lA.ip . 
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B About the uniformization of Q{x,y) = (see |4], chapter 3) 



When the associated Riemann surface is of genus 1, the algebraic curve Q{x, y) = admits 
a uniformization given in terms of the Weierstrass p function with periods iJi,UJ2 and its 
derivatives. Indeed, setting 

Q{x,y) = a{x)y^ + b{x)y + c{x), 
D{x) = b'^{x) — 4:a{x)c{x) = d^x'^ + d^x^ + d2x'^ + dix + do, 
u = 2a{x)y + b{x), 

the following formulae hold (see |i4], lemma 3.3.1). 

(i) If (i4 ^ (4 finite branch points xi, . . . X4) then D'{xi) > and 

D'{xi) 



x{oj) = X4 + 



U[LiJ] 



p{u)-iD"{x,y 

D'{x^)p'{u) 

v2 • 



2{p{u:)-\D"{x,)Y 
(ii) If (i4 = (3 finite branch points xi,X2,X3 and X4 = 00) then 



x{uj) 



d2 

3 



^3 

2d3 ■ 
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